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(Nl . ABSTRACT 



We have used the EFAR sample of galaxies to investigate the light distributions of early- 
type galaxies. We decompose the 2D light distribution of the galaxies in a flattened spheroidal 

■ component with a Sersic ( 1968) radial light profile and an inclined disc component with an 

exponential light profile. We show that if we assume that all galaxies can have a spheroid and 

disc component, that then the brightest, bulge dominated elliptical galaxies have a fairly broad 

distribution in the Sersic profile shape parameter rig, with a median of about 3.7 and a sigma 

JL | of ~0.9. Other galaxies have smaller ns values. This means that spheroids are in general less 

(j) ' i 

5_i ■ concentrated than the de Vaucouleurs R * -law profile, which has rts=4. 

\ While the result of our light decomposition is robust, we cannot prove without kinematic 

C3 • information that these components are spheroids and discs, in the usual sense of pressure- and 

rotation-supported stellar systems. However, we show that the distribution of disc inclination 
angles is consistent with random orientation if we take our selection effects into account. If 

■ we assume that the detected spheroids and discs are indeed separate components, we can draw 
^ ' the following conclusions: 1) the spheroid and disc scale sizes are correlated; 2) bulge-to-total 

luminosity ratios, bulge effective radii, and bulge rig values are all positively correlated; 3) 
the bivariate space density distribution of elliptical galaxies in the (luminosity, scale size)- 
plane is well described by a Schechter luminosity function in the luminosity dimension and 
a log-normal scale-size distribution at a given luminosity; 4) at the brightest luminosities, 
the scale size distribution of elliptical galaxies is similar to those of bright spiral galaxies, but 
extending to brighter magnitudes; at fainter luminosities the scale size distribution of elliptical 
galaxies peaks at distinctly smaller sizes than the size distribution of spiral galaxies; and 5) 
bulge components of early-type galaxies are typically a factor 1.5 to 2.5 smaller than the disks 
of spiral galaxies with a slight luminosity dependence, while disc components of early-type 
galaxies are typically twice as large as the discs of spiral galaxies at all luminosities. 

Key words: galaxies: elliptical and lenticular, cD - galaxies: fundamental parameters - galax- 
ies: statistics - galaxies: structure 



1 INTRODUCTION 

Studying the light distribution of galaxies in order to determine 
and compare their structural properties has a lo ng tradition (e.g., 
Ide Vaucouleursl fl948l : lFisi]|l964t lFreemanlll97(l) . Traditionally, 
the light distributions of elliptical galaxies have been fitted by 
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the R 4 -law profile of Ide Vaucouleursl ll948t) . while spiral galax- 
ies have been fitted with exponential light profiles for their discs 
and often, d ue to their similarities to early-type gal axies, Ri-law 
bulges (e . g . . Ide Vaucouleursll 1 95 9t IKormendyll 1 9771) . In the recent 
decade we have seen some shift from this paradigm. A gener- 
alization of the R 3 and exponential light profile laws described 
first by <Sersidll968l) . resulting in an Rn law, has more often 
been fitted to early-type galaxies to account for observed devi- 
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ations from the Ra profiles (e. g., ICaon. Capaccioli. & D'O nofrio 
1 19931 iGraham & Collesdfl997l) . Based on this development and 
based on theoretical predictions for some forms of bulge formation 
(e.g.. lCombes et al]|l990l : IPfenniger & Normanl fl990). more gen- 
eral forms of light profiles have been fitted to the b ulge light of spi- 
ral galaxies as we ll jAndredakis, Peletier, & Balcellsll995l : lde Jond 
1996a; lMac Arthur. Courteau. & Holtzmanl2003h . With the realiza- 
tion that man y earl y-type galaxies have disc-like c omponen ts (e.g ., 
iBender et al]|l989t iRix & Whitelll99rA Lforgensen & Franxlll994h . 
bulge/disc decompositions of elliptical galaxies have been per- 
formed, but in gene ral with bulges and exponential discs (e.g., 
ISaglia et alJll997bh except for a few notable exceptions on small 
samp l es of galaxi es to'O nofricl200ll : lBalcells. Graham & PeletieJ 
120041 iGutierrez et alJ l2004). Here we perform 2D bulge/disc de- 
compositions on the large EFAR sample of early-type galaxies 
I We g n~et al Jl 1 99(1 hereafter Paper I) using an R~ spheroid and 
exponential disc light profiles, determine relations between their 
structural parameters and their frequency of occurrence, and com- 
pare those to the properties of disc dominated galaxies. 

In an earlier paper in this series <Sagliaetalll997d. Paper III 
hereafter) we presented the luminosity profiles for the 776 galaxies 
observed in the EFAR project to which we fitted seeing convolved 
i?^-law bulge plus exponential disc models. In Paper III, 31% of 
the sample proved to be spiral or barred spiral galaxies and 69% of 
the sample were classified as early-type galaxies. Of those early- 
types 18% were well fit by an R 1 ^ 4 alone, the great majority (70%) 
were very much better fit by the combination of an exponential 
disc plus R± bulge than by a bulge alone. The remaining 12% of 
the early-type galaxies were classified as cD, without reference to 
their profile type, as they were the brightest in their cluster and had 
half-total-light radii larger than lOkpc. 

In Saglia et al. 1 1997b, hereafter Paper IV) we showed that 
the combination of an R% -law spheroid plus an exponential disc is 
degenerate to an R~ -law spheroid, unless one has data for a very 
large range in surface brightness and radius. Given that only 12% 
of the EFAR galaxies were best fitted with an R& law alone, and 
that the remainder of the galaxies were better fitted with an R 3 
bulge and an exponential disc, one might wonder about the validity 
of the Ri law in general. To address this issue, we have refitted all 
EFAR galaxies again with an R~ -law spheroid and an exponential 
disc. If early-type galaxies are dominated by R~ spheroids, most 
fits should converge to largely spheroid dominated systems with 
small discs. If the i?^-law spheroids plus disc picture is correct, 
most fits should converge towards spheroids with Sersic parameter 
n close to 4. Obviously, the truth could be somewhere in between 
with spheroids with Sersic parameter n covering a large range and 
with a large range in bulge-to-disc ratios. 

We have to realize that even when we find that most early- 
type systems have light components that are well described by ex- 
ponential light profiles, this in itself will never prove the existence 
of exponential discs. First of all, there is no reason why the expo- 
nential discs observed in most spiral galaxies should be present in 
a similar way in bulge dominated systems, especially realizing that 
exponential discs in spiral galaxies are still poorly understood. Sec- 
ondly, even if light is present in exponential distribution, it might 
not be in a disc configuration flattened by rotation as discs of spiral 
galaxies are normally understood. One can probably make many 
stable configurations in bulge-like potentials that have small expo- 
nential like stellar distributions on top of them, without the need of 
any net rotation. The best confirmation of exponential light distri- 



butions in early type galaxies being rotationally flattened discs as in 
spiral galaxies has to come from measurements of their kinematics. 

Even when kinematics is availa ble, it still might b e difficult 
to reach an unambiguous result. lScorza & Bended <P995l) made dy- 
namical models for a small sample of early-type galaxies based on 
their photometric bulge/disc decompositions and showed that the 
observed asymmetries in the long-slit stellar absorption line pro- 
files were consistent with the dynamical model predictions. How- 
ever, with their data and modelling technique they could not prove 
this was an unique kinematic interpretation. Better kinematic con- 
straints on discy components can be obtained with full 2D velocity 
profile measure ments, such a s ob tained with the SAURON project 
iFalcon-Barroso et alJl2004l lEmsellem et aljEo04l) . However, so 
far these studies have been limited to the high surface brightness 
central regions of galaxies. 

In recent years a number of studies were performed where 
the galaxy samples were large enough and the selection bias un- 
derstood well enough that bivariate space density distributions of 
structural parameters could be calculated. ICross et alj 1200 lh pre- 
sented the bivariate (lumino sity, sur face brightness distribu tion) for 
the 2dFGRS IColless et alll2001bl) . iBlanton et all <2003l) showed 
many different bivariate space density distributions derived from 
the Sloan Digital Sky Survey (SPSS: lYork et alj|2000h . and pre- 
sented some distributions in bins of colour and concentration in- 
dex to split their sample in different galaxy types. iKauffmann et alJ 
(2003) derived mass-to-light ratios for the SPSS galaxies and pre- 
sented stellar mass bivariate distributions. Ishen et a l] 120031) also 
investigated the bivariate distribution of structural parameters of 
different subsets of SPSS galaxies. However, even though calcu- 
lated, they did not include the absolute space densities in their 
analysis. While based on large sample to provide good statistics, 
in none of these studies bulge/disc decompositions were performed 
to derive space density distributions for the individual components. 
Such an analysis was perfor med on the discs of a s ample of about 
1000 Sb-Sdm galaxies bv Ide Jong & Lace^ l2000h . However, no 
such analysis on early-type galaxies has appeared so far. 

A fundamental question for galaxy formation is the origin of 
the Hubble sequence. Why do some galaxies have discs and others 
not? If in fact almost all galaxies have a disc component but of 
varying mass, luminosity, and size, the question of the origin of the 
Hubble sequence takes on a different perspective; in semi-analytic 
models of galaxy formation early-type galaxies go through phases 
where they have substantial discs which are disrupted and re-built 
during the cycle of hierarchical merging. Petermining galaxy disc 
and bulge parameters and their bivariate distributions will help us to 
test and constrain such hierarchical semi-analytic models of galaxy 
formation. 

In this paper we will use a 2D fitting algorithm to measure 
the bulge and disc properties in the EFAR sample of early-type 
galaxies. We estimate the frequency, luminosity and size of discs 
and bulges and explore the properties of the discs in the context of 
SO and spiral galaxy discs. 

In the next section we describe our method of 2P model fit- 
ting and parameter error estimation and compare our new results 
with the results of Paper III. In section 3 we characterize the prop- 
erties of the discs and their relation to the host bulge component. 
Our conclusions and the implications of this work for models of 
galaxy evolution are summarized in section 4. We use a cosmology 
of fi m =0.3, f2A=0.7, and H =1Q km s" 1 Mpc" 1 throughout. 
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2 GIM2D: 2D GALAXY MODELLING 

We decided to use full 2D modelling of the galaxy images, as this 
provides extra bulge/disc separation constraints due to the differ- 
ent flattening of bulges and discs, sometimes combined with the 
constraints resulting from the different position angles of the two 
components. Performing the modelling in 2D also allows proper 
weighting of data points reflecting uncertainties due to photon 
statistics, flat fielding uncertainties, sky level uncertainties and 
model imperfections (i.e. real galaxies are not as smooth as our 
models). 

To extract the best possible bulge and disc parameters and to 
estimate their errors we performed the following steps: 

Many of the EFAR galaxies have several repeat observations, 
so for the 775 objects more than 2500 images are available. We only 
used the 7?-band images, because the B-band imaging of the EFAR 
sample is incomplete. We determined the seeing on all EFAR im- 
ages by finding isolated star-like objects in each of them, then mea- 
sure the FWHM of each object by fitting a Gaussian profile, and 
taking the error weighted average of the eight objects in each im- 
age with the smallest errors after removing the objects with smallest 
and largest FWHM. At this stage we selected for each galaxy the 
image with the best seeing. All galaxies with too few good fore- 
ground stars (< 3) to make a reliable seeing measurement were 
removed from the sample (23 galaxies). 

SExtractor iBertin & Amoutsl 19961) was used to create a mask 
of all objects on the images, splitting the images of overlapping ob- 
jects in separate components. This mask is used by the fitting rou- 
tine to exclude all pixels with objects detected except for the pixels 
with the galaxy of interest and the sky pixels. A very low level SEx- 
tractor mask was also used to determine the sky level of the image, 
using a kappa-sigma clipped median of the pixels deemed free of 
objects, and this sky level was kept fixed during the fitting of the 
galaxy. 

The 2D modelling of the galaxy images was performed by 
the publicly available GIM2D fitting package (Version 2.2.1), de- 
scribed in detail in Simard et al. (2002). The spheroidal component 
of each galaxy was fitted with a Sersic law: 

SeM^e-^^] 1 /" 13 ^ (1) 

where the effective radius (-R c ,b) encloses half the total luminosity 
and E e is the surface brightness at this radius. The parameter c is 
set to 1.9992n_g -0.3271 so that i? e ,B remains the half total light ra- 
dius of the bulge. Better approximations of c were published after 
we d id our calculations (e.g., MacArthur. Courteau. & Holtzman 
2003), but these improvements are only relevant for tib values 
much less than 1, which are rare in our sample. The bulge com- 
ponent was allowed a flattening with ellipticity e = 1 — b/a, with 
a, b semi-major and minor axis. The i? c ,B provided by GIM2D is 
the bulge half-total-light radius of the major axis, which we mul- 
tiplied by VI — e to be consistent with our previous work, which 
uses the effective radius on a circular aperture. 

The disc was modelled with an exponential disc: 

E D (r) = E e-^ h , (2) 

which was assumed to be infinitely thin, transparent and have an 
inclination i. The two galaxy components had a common center 
which was fitted as well, and an independent flux, here expressed 
in a bulge-to-total light ratio (B/T). Each component has also 
an independent position angle. The summed galaxy components 
were convolved with a Gaussian with a FWHM equal to the de- 
termined seeing, using an FFT algorithm. GIM2D fits this large 



1 1 -dimensional param eter space of the two components with the 
Metropolis algorithm ^Metro polis et alj[l953l) . The Metropolis al- 
gorithm is very CPU intensive, but it tends to be more robust than 
methods based on gradient searches. The Metropolis algorithm is 
ideal for cases as the current where many fits have to be performed 
unsupervised. The CPU intensiveness was solved by using 16 fast 
Linux Dual Processor PCs in a cluster. 

The fitted scale size parameters in pixels and intensities in 
counts were converted to arc-seconds and magnitudes using cali- 
brations described in Paper III. To set the nomenclature, the galax- 
ies are parametrized in the remainder of the paper as follows: for 
the total galaxy we have luminosity M to t, effective (half-total- 
light) radius R c and effective surface brightness within this radius 
(SB C ); for the bulge luminosity Mb, effective radius R c ,b and ef- 
fective surface brightness (S-B ,b); for the disc luminosity Mo, 
scale length h and central surface brightness /io- 

After all galaxies had been fitted, galaxies with high \ 2 resid- 
uals and/or suspicious model subtracted residual images were iden- 
tified. In many of these cases SExtractor had not detected objects 
near the nucleus of the target galaxy or done an insufficient splitting 
of nearby objects, resulting in incorrect masks and causing poor 
fits. For most of these objects we were able to modify the masks by 
hand to give satisfactory result, but a number of objects with bright 
foreground stars, many bright companions or the occasional galax- 
ies with double nuclei were removed from the sample at this stage 
(20 galaxies). 

We used the model subtracted residual images to identify 
galaxies with spiral structure. Even though not necessarily spiral 
galaxies based on their B /T ratio, these galaxies are marked with 
special symbols in the diagrams. We will use the term "non-spiral 
galaxies" in the remainder of this paper for those galaxies that had 
no clear spiral structure in the residual images (454 of the 558 
galaxies used in the final selection). The brightest galaxy of each 
of the clusters with at least 3 mem bers as identified in Paper VII 
of this series IColless et al. 2001a) is also marked with a special 
symbol in the diagrams. 

The GIM2D fit package provides 68% confidence limits on 
the fitted parameters, based on the topology of the parameter space 
that is build up during the fit process. These errors take only the 
uncertainties in the fitted parameters into account, not uncertainties 
in fixed parameters like sky level and seeing. As errors in these pa- 
rameters ar e often the main error in the determined bulge and disc 
parameters I de Jong 1996a), we repeated all fits with the maximum 
error expected in our determination of sky and seeing values. 

We estimated the typical error in our sky estimates by dividing 
the image in 4 and measuring the sky in each of the quadrants sepa- 
rately using the same procedure as for the total image (using object 
masking and kappa-sigma clipping). The RMS in the 4 sky aver- 
ages provides an estimate in the sky level uncertainty on a galaxy 
size scale. This is probably a slight overestimate, as most galaxies 
are near the center of the frame, so near the average sky of the total 
frame, while the sky RMS of the quadrants will be more affected by 
sky level gradients across the frame. Most images used have small 
sky level uncertainties: 80% of the galaxies have quadrants errors 
less than 0.5%, only 8% of the galaxies have quadrant RMS values 
larger than 1% of the sky, mainly due to sky gradients. We repeated 
all GIM2D fits with 1 % sky added and subtracted, and the resulting 
fit parameters should be seen as the maximum errors in the param- 
eters that can be expected to be caused by sky uncertainty. 

For each galaxy we repeated the fit with seeing increased and 
decreased by 2% to determine the errors in the fitted parameters due 
to seeing measurement uncertainties. Judging from the RMS in the 
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seeing measurements of different stars in the different frames, less 
than 20% of the galaxies have uncertainties in their seeing larger 
than 2%. Again, these errors should be seen as the maximum error 
in the fitted parameters due to seeing measurement uncertainties. 

Where appropriate, we show typical errorbars on our graphs, 
which are the median errors for all galaxies in the determined pa- 
rameters, where the parameter error for an individual galaxy is the 
maximum error due to fit parameter space, sky uncertainty, or see- 
ing uncertainty. We chose the maximum error here, because it is 
hard to assign relative weights to the different error sources. To 
assign proper weights would require a full set of Monte Carlo sim- 
ulations, which is prohibited given the large computation times in- 
volved. 

We excluded all galaxies with no measured redshift from our 
analysis, mainly spiral galaxies which were excluded from the 
spectroscopy phase of the EFAR project as being not interesting 
for the main aim of the project. The final sample of galaxies used 
contains 558 galaxies. 
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2.1 Testing the 2D fit routine 

We tested the GIM2D fitting program by performing fits on 500 
artificial galaxy images with characteristics similar to our real data. 
The images were 500x500 pixels and the pixel size was 0.65 arcsec. 
Skylevel was set at 20.5 mag arcsec" 2 . Galaxy images were created 
in a Monte Carlo fashion with total magnitudes sampled linearly 
between 12 and 16 mag, with B/T ratios sampled linearly between 
0.05 and 1, and having bulges with 18 < (SB C , B ) < 22 and 1 < 
nfl ^ 6.5, and discs with 20 ^ /Uo S= 23. CCD readnoise and gain 
were set comparable to most of our observations and Poison noise 
was added appropriate for the detected flux level. The galaxies were 
convolved with a Gaussian seeing with FWHM of 1.8 arcsec, which 
is at the high end of the typical seeing of the used images (90% with 
seeing less than 1.8 arcsec). 

A comparison of the fitted GIM2D parameters with the in- 
put parameters is presented in the first five columns of Fig.Q All 
parameters are reasonably well reproduced, but some systematic 
trends can be seen. As could be expected, the largest errors in the 
bulge parameters (i? c ,B, «b) occur at small B/T values, while the 
largest disc errors can be found at large B/T values and small h 
values. 

Obviously, for real galaxies we do not know the "input" pa- 
rameter values, but examining the errors in the fitted parameters as 
function of the fitted output parameters, shows parameter ranges 
we should not trust. We therefore exclude here from the artificial 
data and later from the real data all fitted R c ,b values for galaxies 
with B/T <0.2 or i? c , B < 1-0 FWHM seeing. Likewise, we exclude 
all fitted h values for fits where B/T>0.95 or h<l.5 FWHM see- 
ing. Applying these cuts results in RMS errors between input and 
fitted parameters of order 5% in M tot , 10% in B/T, 30% in R C , B , 
20% in n B , and 15% in h. 

It should be noted that some of the errors are strongly corre- 
lated, as shown in the last two columns of Fig.Q Especially the 
bulge parameters i? e ,B and ns are strongly correlated, with some 
weaker correlations between these bulge parameters and Mtot and 
B/T. 

The ub values show a systematic trend with the input ns in 
Fig.Q This is for a considerable fraction the result of the limits 
imposed on the fitted ub values, which were constrained to 1 ^ 
ub ^ 6.5. This means that artificial galaxies with ub values close 
to 1 will have fitted ns values preferentially scattered to larger ns 
fit values, and likewise for galaxies with input ns values close to 



Figure 2. The integrated magnitude versus effective radius difference com- 
paring the Paper III to the current GIM2D determinations. The effective 
radius difference is defined as AiJ c = 2(R C 1D - R C GIM2D ) / (R C 1D + 
R e GI M2D ). Open square denote brightest galaxies in EFAR clusters, star 
symbols denote galaxies with clear spiral structure in their model subtracted 
residual images. The errorbars indicate the median error for all galaxies in 
the parameters, where the parameter error for an individual galaxy is the 
maximum error due to 2D fit parameter space, sky subtraction error, and 
seeing error. 



6.5. A similar effect can be seen in the AB/T versus B/T plot. 
To minimize this effect on real galaxies we increased the limits to 

1 sC n B < 8. 

At small ns values a relatively large number of galaxies show 
"catastrophic" errors. These are cases of mistaken identity; the 
bulge and disc shapes are so similar that they are swapped. Without 
prior assumptions about relative scale sizes or surface brightnesses 
this is unavoidable. 



2.2 Comparing profile fitting and 2D fitting 

In Paper III we derived bulge and disc parameters of the EFAR 
galaxies using ID luminosity profiles. Circularly averaged lumi- 
nosity profiles were fitted by the combination of an R 3 and an 
exponential luminosity light distribution, convolved with seeing, 
which was also fitted. Multiple exposures of the same object were 
optimally combined and an optional sky-fitting procedure was used 
when necessary to correct for sky subtraction errors. 

We will now compare the Paper III ID parameters with the 
current GIM2D parameters to study the differences between ID 
and 2D fitting and R& and _R«-law bulges. We removed the fit 
quality Q=3 fits from the Paper III results, as those were considered 
of poor quality. 

Figure|2| shows the difference in integrated magnitude ver- 
sus the effective radius difference (defined as AR C = 2(R C 1D — 
R a GIM2D )/(R c 1D + R GIM2D )) when comparing the ID and 
the GIM2D results. The GIM2D magnitudes are 0.06 mag fainter 
in the median and agree to within 0.17 mag RMS. The effective 
radii agree to within 28% RMS, without taking the strong correla- 
tion into account. The correlated differences are for a considerable 
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Figure 1. Results of GIM2D fits to the Monte Carlo simulated artificial galaxy images. The first five columns on the x-axis show the input parameters of the 
galaxy images, the y-axis and the last two columns on the x-axis show the differences between fitted and input parameters, defined as: AMtot = Af£t ~ -^tof 
AB/T = B/T &t - B/T in , An B = n|* - rag, Ai?, c , B = log(R^ B /R^ B ), and Ah = log(/i flt /h in ). 



fraction due to the change in bulge parameter shape allowed by the 
R~ law. Even though the flux in the inner region of the galaxies 
is well determined, the total flux (and hence half light effective ra- 
dius) is determined by extrapolating the model to infinity. This is 
demonstrated more clearly in the next figure. 

Figure|5| shows the old versus new bulge effective radius dif- 
ference as function of Sersic bulge shape parameter n b ■ When the 
2D fitted ub is getting smaller than the n B —4 value, light is 
extending less far and R c ,b is getting smaller in the R~ law case 
than in the i?3 law case. The opposite happens when ns > 4. 



Taking this effect into account, the bulge parameters of the galax- 
ies agree to within 50% RMS of each other, a larger error than for 
the total galaxy parameters as might be expected. 

Figure |4] compares the old and new disc scale lengths. Some 
of the extreme outliers are due to the swapping of interpretation of 
which component is disc and which is bulge. What used to be a 
small disc in a large i?3-law bulge has become a small _R~-law 
bulge with in general low n_g and sometimes vice versa. The disc 
scale lengths agree with each other to almost 35% RMS, which sug- 
gests that many early-type galaxies seem to have a light component 
that is robustly fitted by an exponential light profile, to some extent 
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Figure 3. The bulge effective radius difference between ID Paper III deter- 
minations and the current GIM2D determinations as function of bulge shape 
rig. The effective radius difference is defined as AR C b = 2(-R c ,b 113 — 
Rc,B GIM2D )/(Rc.B 1D + R c ,b GIM2D )- Symbols' and errorbars as in 
Fig.H 



Figure 4. The disc scale length difference between ID Paper III deter- 
minations and the current GIM2D determinations as function of bulge 
shape rig. The effective radius difference is defined as Ah = 2(hj D — 
h GIM2Dy, h lD + h GIM2Dy wherf , h lD = h W / cos(i GIM2Dy 
Symbols and errorbars as in Figfil 



independent of the shape of the spheroid component. This does not 
prove that these are real discs: other radial disc light profiles may be 
more appropriate (see e.g. lScorza & Bendedll995|) and kinematics 
is required to confirm a disc configuration for an individual galaxy. 

In this section we showed that while the ID and the 2D 
R n fit result in systematically different structural parameters, these 
deviations are correlated with the Sersic shape parameter. In deter- 
mining uncertainties in the parameters we have to take these sys- 
tematic chances into account. Therefore, under the assumption that 
the combination of an R ~ law spheroid and exponential disc is 
a reasonable good description of the intrinsic light distribution of 
a galaxy and by giving some more weight to the artificial images 
test of Section 12.11 we conclude that bulge and disc parameters can 
be determined with a typical 1 sigma accuracy of about 30%. How- 
ever, one should keep in mind that some extreme outliers with much 
larger errors will be common in large samples. In the remainder of 
the paper we will use the GIM2D determined parameters to investi- 
gate correlations between bulge and disc parameters and the space 
density of early-type galaxies as function of these parameters. The 
fitted GIM2D parameters are available in electronic format from the 
Centre de Donnes astronomiques de Strasbourg (CDS) web site. 



3 PROPERTIES OF DISCS AND BULGES 

3.1 Sample selection and selection effects 

The sele ction of the EFAR galaxy sample has been described in 
detail in [Wegner et alJ ll996l) . and we will repeat here the most 
important aspects relevant for the current study. All galaxies were 
selected with the intention to be cluster members, with the clusters 
lying in two selected regions in the sky which are rich in galaxy 
clusters. The target clusters had redshifts in the range from 6000 
to 15000 km s" 1 . Photographic enlargements of POSS and SERC 



Sky Survey plates of the target clusters were inspected for galaxies 
with elliptical/SO galaxy like morphology and a minimum diameter 
of about 16 arcsec. Morphological selection was done to err on the 
safe side, including many spiral galaxies which were identified as 
such later with subsequent CCD imaging. The sample is therefore 
intended to be complete in large, early-type galaxies (elliptical/SO), 
but is not complete in any other galaxy type. In most of the analysis 
presented here, we will only use the galaxies from the EFAR tar- 
get clusters, excluding the galaxies from the comparison sample in 
Coma, Virgo and the field, as they suffer from completely different 
selection effects. 

To maximize the number of early-type galaxies, galaxies were 
selected with rather round isophotes and the sample is therefore 
rather face-on. The sample has only eight galaxies with ellipticity 
e>0.4 at the outer isophotes. The distribution of intrinsic shapes 
of early-type galaxies is poorly known, and we have therefore not 
attempted to correct the bulge parameters to face-on values using 
the bulge isophote shapes. However, the disc parameters shown in 
the diagrams are face-on values, assuming fully transparent discs. 

While we can not prove that our fitted disc light distributions 
are really discs, we can use the minor over major axis ratio (bla) 
distribution to make a statistical test. For a sample of infinitely thin 
discs at random viewing angle the distribution of bla is expected to 
flat. We show this distribution for all non-spiral galaxies in Fig.|5| 
(line). The distribution is rather flat for b/a > 0.6, but steeply de- 
clines for smaller values of bla. This is the result of our selection 
bias, with very few galaxies having e>0.4 (i.e. 6/a<0.6). Assum- 
ing that discs are in general aligned with the spheroid they are sit- 
ting in, we expect the bla distribution to be peaked toward face-on. 
This effect will be less for small discs in large bulges and indeed 
the distribution becomes relatively more spread out when looking at 
all non-spiral galaxies with B /T>Q.l in Fig.[5](solid histogram). 
While not conclusive, this test shows that our fitted disc compo- 
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Figure 5. The distribution of minor over major axis ratios (bla) for all 
non-spiral galaxies (line histogram) and non- spiral galaxies with B/T>0.1 
(solid histogram). 

nents have orientation angles consistent with the random distribu- 
tion expected. 

To fully quantify the selection limits of the sample as function 
of the structural parameters of the galaxies, we need not only to 
know the angular diameter limits, but also the surface brightness 
at selection diameter. The average surface brightness at the selec- 
tion diameter was 22.05 i?-mag arcsec" , with a RMS of 0.5 R- 
mag arcsec" 2 , as derived from the ID profiles of Paper III. This 
value is used, in combination with the typical log(Z?w)>1.2 diam- 
eter limit (Paper I, see also Eq.|3J and the minimum 6000 km s -1 
distance limit, to calculate the selection limits indicated in our dia- 
grams where possible. 

3.2 Galactic extinction and cosmological corrections 

The observed parameters were corrected for Ga lactic extinc- 
tion according to the precepts of iBursteinl 1200 31). which aver- 
ages the estimates of [ Schlegel. Finkbeine r. & Davisl 1 19981) and 
lBurstein&HeilesNl98a with appropriate offsets. No attempt was 
made to correct the parameters for internal extinction. We used 
the cluster redshifts as listed in Paper VII to calculate distances 
to individual galaxies. Cosmological corrections were made using 
fi m =0.3, S7a=0.7, and Ho=70 km s _1 Mpc -1 , and K-corrections 
were made as described in Paper III. All these corrections are typ- 
ically small for these galaxies and making different reasonable as- 
sumptions will not affect our results at all. 

3.3 Structural parameter correlations 

We will now investigate correlations in structural parameters of 
bulge and disc. Obviously, many combinations of bulge and disc 
parameters can be made, and we will show here only those deemed 
most interesting. 

In Fig.|6|we show bulge effective radius versus bulge effective 
surface brightness for our sample. There seems to be a reasonable 
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Figure 6. Bulge effective radius versus bulge effective surface brightness. 
The solid line shows the selection limit for pure R * profile galaxies at 6000 
fans -1 . The dashed line indicates Ri profile galaxies of -25 absolute R- 
mag. Bulges with equal luminosities will lie on lines parallel to this line. 
Symbols and erorbars as in Fig.l2l 

tight correlation between these two parameters for the non-spiral 
galaxies, but we have to be wary of selection effects. The solid 
line in this diagram shows the selection cutoff for galaxies with an 
Ri profile dominating at the selection radius, calculated using the 
average surface brightness at the selection radius. The rather sharp 
cutoff in the high surface brightness, long scale size region of the 
diagram is therefore real; there are no selection effects acting in 
this part of the diagram. This upper limit is related to the zone of 
exclusion in K-space discussed in detail by Bur stein et aljil997t) . 

At the low surface brightness end of the diagram we have 
to be more careful: most of the galaxy distribution seems to 
lie clearly separated to the right of the selection limit line, but 
the selection line is for a 6000 km s -1 cluster, and moving to 
the right for the higher re dshift clusters. The "visibility volume" 
toisnev & PhirliDPsll983h for these lower surface brightness galax- 
ies is definitely much smaller than for high surface brightness 
galaxies, certainly taking further into account the measurement un- 
certainty in the selection diameter (Paper I). We will make full "vis- 
ibility volume" corrections in Section l3~4l to calculate the volume 
density of early- type galaxies as function of structural parameters. 

There are quite a number of galaxies, whose bulges do 
not follow the main trend seen for the most non-spiral galax- 
ies in Fig.[6] These bulges lie often to the left of the selec- 
tion line, and these galaxies are included in the sample be- 
cause their discs make a significant contribution at the selec- 
tion diameter or because their profile does not follow an R? 
law. Many of these bulges, especially in some of the spiral 
galaxies, may have the more exponential-like bulges typical of 
late-type spiral galaxies | Andredakis, Pele tier. & Balcellsl 1 19951 
ICourteau. de Jong. & Broeilslll996llGraharnll200ll) instead of the 
more R* like bulges found in early-type galaxies. 

In Fig.Qwe show the disc scale length versus disc central sur- 
face brightness for all fitted galaxies. The dotted line shows the 
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Figure 7. Disc scale length versus disc central surface brightness. The 
dotted line shows the selection limit for exponential discs at cz = 
6000 km s — 1 . The dashed line indicates exponential discs of -25 absolute 
ij-mag. Symbols and errorbars as in Fig.l2l 



Figure 8. Disc scale length versus bulge effective radius. To guide the eye, 
the line of equality is shown as a dashed line. Symbols and erorbars as in 
FigJl 



cz — 6000 km s _ selection limit for galaxies which have an ex- 
ponential light distribution dominating at the selection radius (i.e. 
spiral galaxies). Elliptical galaxies will not adhere to this selection 
limit, as they were selected on basis of their R~ light distribution. 
Some interesting trends can be seen in this diagram, but we have 
to be wary of more hidden selection effects. At the very short scale 
lengths we are close to our scale length cutoff, as 1 kpc is of order 
2.5 arcsec for cz — 6000 km s _1 and even less for higher redshifts. 
The observed limit at the low surface brightness, small scale length 
end is more subtle and probably arises from the minimum contrast 
needed to see a low surface brightness disc on top of a bright bulge 
distribution. Using the relation between h and R c ,b shown in Fig.|8| 
discussed in the next paragraph and the relation between J? c ,b and 
(SB c ,b) shown in Fig.|6| we find that we probe discs with central 
surface brightnesses at most about 2 magnitudes fainter than the 
(SB C: b) of the bulge. It is not unreasonable to suspect that the low 
surface brightness cutoff is caused by lack of disc/bulge contrast. 
The only cutoff observed in this diagram that is definitely real is at 
large scale length, high surface brightness, as no selection effects 
operate in this part of the diagram. This cutoff is similar to the one 
observed for spiral galaxies ide Jons 1996b). 

In Fig.|8|we show the correlation between bulge and disc scale 
sizes. No obvious selection effects operate in this diagram, except 
for a spiral galaxy having a large disc scale length and a small bulge 
would probably not have entered the sample. This is not the case 
for elliptical galaxies and therefore we have to make sure that other 
effects are not causing this correlation, before we can accept it as 
real. Our Monte-Carlo simulations showed that our GIM2D fitting 
routine is capable to recover models of the full parameter space 
presented here, so it is unlikely that it is the result of a fitting rou- 
tine artifact. Indeed, we do see a small number of galaxies with 
R e ,B/h > 1 in Fig.|8| and so the routine is capable of covering this 
parameter space in real data. Still, given the limits on the parameter 
space for a given galaxy (not larger than image size, not smaller 



than seeing), it is not completely surprising that the bulge and disc 
scale parameters are similar in size. 

The R c ,B/h > 1 galaxies can be divided in three groups. A 
number of these galaxies are clearly disturbed by dust-lanes or are 
poorly fitted due to nearby bright galaxies or stars. We may have 
a number of cases of mistaken identity, where bulge and disc are 
swapped. This is especially the case when is low and the bulge 
has an almost exponential light profile. Swapping the bulge and disc 
parameters on these galaxies would bring them in line with the gen- 
eral trend. Finally there are a very few galaxies which seem to have 
a genuine small inner light component for their given bulge size, but 
these should not be mistaken for the nuclear discs found in HST im- 
ages of early-type gal axies, as these a re a factor 5-20 smaller yet 
again iScorza & van den Bosch 1998). In a number of cases they 
might be identified with the inner break radius observed in HST im- 
ages of earl y-type galaxies jFaber et alll997ilTruiillo et alj2004h . 
as shown bv lD'Onofrid 1200 ll) . 

In Fig.|5|we show in three sets of diagrams the bulge, disc, 
and combined bulge and disc parameters as function of B/T ratio 
and Sersic ub values. Before we can draw any conclusions on any 
apparent correlation, we have to be aware of the many selection 
effects playing in these diagrams. Some are explicit, like the limits 
we placed on some parameters based on B/T ratio to reduce fit 
errors. Others are more hidden, like the exclusion of bright discs, 
as their more visible spiral structure would have precluded them 
from entering the sample in the first place. 

Concentrating first on the B/T ratios, we see that more bulge 
dominated systems have larger bulge scale sizes, but lower effec- 
tive surface intensities, which combines in a rather weak B/T ratio 
trend with total bulge luminosity. On the other hand, looking at the 
disc parameters in the middle set of diagrams, we see a weak trend 
of disc scale sizes with B/T ratio, but at each B/T ratio there 
is a wide range in disc scale sizes. More bulge dominated systems 
have discs of lower surface intensities. The combined trends of disc 
scale size and surface brightness result in a strong correlation be- 
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Figure 9. Bulge-to-total-light ratio and Sersic profile shape parameter nj as function of the bulge parameters (top), disc parameters (middle), and bulge- 
disc combined parameters (bottom). Symbols and errorbars as in Fig.|2] The dashed line in the h/R a ^ versus B/T diagram indicates the expected relation 
between these parameters under the assumption that there is a fixed bulge to disc surface brightness difference with arbitrary zero-point. 
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tween disc luminosity and B /T ratio for our sample dominated by 
bright, early-type galaxies. 

Combining these results for this sample, the change in B/T 
ratio seems to be mainly driven by changes in disc luminosity, not 
bulge luminosity. This is opposite to w hat has been o bserved for 
disc dominated systems as found b y e.g. lTruiillo et alJ l2002t) and 
iBalcells. Graham & Peletierl 120041) . who showed that most of the 
change in B/T ratio in these systems is due to bulge luminosity 
variation. We have to be cautious though. Our selection criteria re- 
sulted in a very limited range in total galaxy luminosities. We are 
only sampling the upper decade in the luminosity function, the 5- 
10 brightest galaxies in a cluster. It could be that lower luminosity 
elliptical galaxies reveal often equal sized discs, changing the ob- 
served trends. 

Finally, we see that the weak opposite trends in the bulge and 
disc scale sizes with B/T ratio combine in a strong correlation 
between /i/7? c ,b and B/T. We should be cautious not to see too 
much in this relationship. The expected relation between h/R e: B 
and B /T under the assumption that there is a fixed bulge to disc 
surface brightness difference and no significant bulge profile shape 
change is indicated by the dashed line. The bulk of the galaxies 
have indeed a limited range in disc to bulge surface brightness dif- 
ference as shown in Fig.|5| It is just hard to see a faint component on 
top of a bright component, unless their light distributions are com- 
pletely different. The variation from tib = 1 to ns — 5 amounts to 
only a factor 2 in bulge luminosity change at fixed bulge to disc sur- 
face brightness ratio (see Graham 2001, his fig. 4) and this change 
is systematic with h/R Ct B as shown in Fig.|5| All in all, the relation 
between B/T and h/R C] B with the observed scatter is expected, 
given the observed scatter in disc to bulge surface brightness con- 
trast and the limited and systematic effect in ns- 

We note that lD'Onofriol 1200 ll) finds a different slope in the 
B/T versus /i/J? c ,b relation. Many of his discs were fitted to the 
small inner components (either discs or cores) of the galaxies in his 
nearby sample. We are not able to resolve these components and, as 
D'Onofrio rightly indicates, it is unlikely that these inner compo- 
nents represent similar structures as the more extended components 
fitted in discy galaxies. We do fit the more extended components 
and indeed find a they have different relation to the main galaxy. 

If we now turn our attention to the relationships between tib 
and the bulge and disc parameters, we note that tib has a behavior 
very similar to the B /T ratio with respect to these bulge and disc 
parameters. There are some more outliers in the ns relations, most 
of which are for ns values larger than 5.5. 

The very similar behavior of ns and B/T ratio in Fig.|9|sug- 
gests that these two parameters must be correlated, which is indeed 
the case as shown in Fig. 1101 As mentioned before, all other things 
being equal a change in ub from 1 to 5 only results only in a factor 
2 change in bulge to disc luminosity, so the change in n_g does not 
drive the variation in B/T ratio on its own. For this set of early- 
type galaxies, a large fraction of the B/T change seems to come 
from the change in bulge effective radius, and it is the correlation 
between n_g and R c ,b that seems to drive the trend between ns 
and B /T ratio. Still, cause and effect are hard to disentangle based 
on this photometric data set alone. We note that these trends can- 
not solely be the result of the correlated errors noted in Fig.Q as 
the slope of the errors is incorrect and the size of the effect is too 
small. 

In Fig.[IT] we show the observed distributions of the ns pa- 
rameter for different subsets of our sample. Galaxies that were fit- 
ted with n b values larger than 6 were outliers in previous diagrams 
and are considered of dubious quality. They are excluded from fur- 
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Figure 10. The Sersic tib bulge shape parameter versus B/T ratio. Sym- 
bols and errorbars as in Fig.l2l 




Figure 11. The distribution of Sersic profile shape parameter ub for our 
sample. The top solid line shows the total sample, the hashed histogram 
shows all galaxies without clear spiral structures and the solid histogram 
show all non-spiral galaxies with B/T > 0.5. 



ther analysis here. The distributions are quite narrowly peaked, with 
1 sigma dispersions ranging from about 1 .05 for the total sample to 
0.88 for bulge dominated non-spiral galaxies. This width is about 
twice the mean of the individual maximal error estimates of the 
ns values, and is therefore significant. The median ns value shifts 
from 3.24 for the total sample, to 3.31 for the non-spiral sample, 
and to 3.66 for the non-spiral, bulge dominated B/T > 0.5 sam- 
ple. For galaxies with B/T > 0.7 the me dian is 3.99. The real 
n b values might be even somewhat lower. IrTalcells et all 120031) 
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showed that for a sample of early type disc galaxies the ns values 
decreased compared to ground based data when a nuclear compo- 
nent revealed in their HST data was included as a seperate compo- 
nent in their fits. 

This result complicates bulge-disc decompositions of higher 
redshift galaxies with typically lower signal-to-noise and resolu- 
tion. On such lower quality data a Sersic-law bulge with expo- 
nential disc does not reliably converge and one is force to use 
a fixed ns value. Given that in this morphologically preselected 
and bright galaxy sample there is a broad distribution of ns inde- 
pendent of luminosity (see Fig.|9j peaking at values below ns=4, 
and given that spiral gal axies tend to have bulges with n_a values 
in the 0.5-3 range ( e.g. jMacArfhur. Courteau. & HoltzmartEooll : 
lBalcellsetaH E()03). it may be more appropriate to use a ns in 
the 1.5-3 range rather than the more customary ns=4 in fixed ns 
fit s. A similar con clusion was reached for early-type disc galaxies 
bv lBalcells et"aT] <2003l) . but is now extended to bright, early-type 
cluster galaxies here. 



3.4 Space density distribution functions 

Probably more important than correlations between structural pa- 
rameters is the true space density of galaxies as function of struc- 
tural parameters, i.e. the distribution of structural parameters cor- 
rected for selection effects. 

When calculating space densities of galaxies from the EFAR 
sample there are several points to note. The sample is not com- 
plete for spiral galaxies, so we can not say anything about space 
densities of those galaxies. The regions selected by EFAR where 
particularly rich in clusters, and therefore are not representative for 
the average space density of galaxies in the local universe. Assum- 
ing that galaxies have similar properties in different clusters, we 
can still use the relative space density of structural parameters, but 
for absolute space densities we will have to re-normalize the EFAR 
regions to the universal cluster density. 

Calculating the space density of galaxies requires good un- 
derstanding of the selection function. The selection of the EFAR 
sample is described in detail in Paper I. The selection of EFAR 
galaxies was a two step process, where first clusters were selected, 
believed to be in the right redshift range and to be rich enough to 
obtain accurate distances. In the next step galaxies were selected in 
each cluster according to a minimum angular diameter crite rion. 

W e will use here a simple Vmax correction method ISchmidtl 
Il968t) to calculate the relative space density of E/SO galaxies 
as function of their structural parameters. The Vmax correction 
method assumes that the average space density of an object with 
certain parameters is inversely proportional to the maximum "vis- 
ibility" volume Vmax of that object, given the selection criteria 
of the sample (for mo re detailed descriptions see lFelterJll976t 
iDisnev & Phillippsll983llde Jong & Lacevl2000l) . The selection of 
galaxies was a two step process and therefore calculating the max- 
imum visibility volume of a galaxy must reproduce this two step 
process. 

First, the Vm ax of each individual galaxy was calculated, by 
integrating the total volume occupied between the minimum and 
maximum selection redshifts of z m i n =6000/c and z max =15000/c, 
with c the speed of light. While doing the volume integral, we have 
to take into account the diameter selection function determined in 
Paper I for each cluster as applied to each galaxy individually. The 
selection function describes the selection diameter (Dw) depen- 
dent probability of sample inclusion, parametrized by: 



P{\og{D w )) = 0.5(1 + erf[log(D w /D%)/5 w ]), 



(3) 



where D\ v is the midpoint and Sw the width of the cutoff in the 
selection function. We furthermore took into account that the con- 
tribution to the volume has to drop to zero when R c becomes less 
than 1.5 arcsec had the galaxy been at that redshift, because we 
excluded such galaxies from our analysis. The integral to calculate 
the total visibility volume is then: 



Vmax = 47Tf2 



D M {z) " (l + z)H(z) 



with 



S(z) 



if 7?, 



Dm ( z ci) 



otherwise 



> 1.5 arcsec 



(5) 



zci the cluster redshift of the galaxy, Qf = 0.139, the fraction of 
the sky covered by our survey, and Dm(z) comoving distance and 
H (z) Hubble constant at redshift z as defined in Hoea 1 1999). 

In the second step, we have to take into account that only ob- 
jects were selected if they belonged to a cluster, with at least three 
early-type galaxies obeying the diameter selection criterion. There- 
fore, the Vmax of a cluster is identical to the V max of the galaxy 
ranked third in diameter belonging to the cluster. The final V ma x of 
each galaxy is the minimum of its individual and cluster V ma x- This 
last step only modifies the Vmax of the top 3 diameter galaxies in 
each cluster, and only if the diameter limit of the cluster is reached 
before the third ranked galaxy is redshifted to 15000 km/s. 

As mentioned before, the EFAR regions were chosen to be 
particular rich in clusters at the target redshift range. We used NED 
to estimate the size of the over density. NED lists 140 original, non 
suppl ement Abell clusters lAbelllfl958t lAbell. Corwin. & Q lowin 
1 19891) for |6| > 30° (fi/ = 0.5) with 6000 < v < 15000 km/s 
versus 87 in the EFAR regions (Qf = 0.139). This translates in 
an over density of about 1.75. However, we observed only 32 of 
the 87 clusters identified by NED, and while not every NED cluster 
should be in our sample (some cluster redshifts are based on only 
one galaxy redshift), it is clear that the over density is probably 
offset by incompleteness. Our sample contains another 41 clusters 
selected by other means, which are not included in the Abell cata- 
log. We have no means to access the completeness of these in gen- 
eral smaller clusters. Given these uncertainties, we have decided 
not to correct for over density and incompleteness, and assume that 
the uncertainty in the absolute space density zero-point is a factor 
of two. We do however correct for our ellipticity incompleteness, 
assuming that elliptical gala xies with e<0.4 con stitute 65% of the 
total E+S0 population l J0rgensen & Franx 1994). 

We will show the factor 2 uncertainty in absolute zero-point 
of the space density distributions in the diagrams presented in the 
remainder of this section. This uncertainty will not affect the con- 
clusions drawn from these diagrams. We will also show that our 
space density distribution and normalization is remarkably consis- 
tent with that of the SDSS, giving confidence that our normalization 
is indeed correct within a factor of 2. 

In Fig.E2l we show the combined bivariate space density of 
all non-spiral galaxies as function of absolute magnitude and effec- 
tive radius. The 95% confidence limit error estimates on the num- 
ber density distribution and the 95% confidence upper limits in the 
regions with no galaxi es in the sample were de rived in a similar 
fashion as described in Ide Jong & Lacevl fcOOOl) . assuming a per- 
fect i?4 law luminosity profile for all galaxies. These error esti- 
mates and upper limits include the effects of distance and diameter 
uncertainties and Poisson statistics. 
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Figure 12. The space density of all non-spiral galaxies as function of total 
magnitude and effective radius of the whole galaxy is shown as the solid line 
histogram. Errorbars on the histogram show the 95% confidence limits due 
to uncertainties in selection diameter and distance and due to Poisson statis- 
tics. The upper limit symbols indicate the 95% confidence limits based on 

i 

the non-detection of such galaxies and the visibility volume of R 4 profiles 
shape galaxies for our selection criteria. The errorbar in the top-left corner 
of the top panel shows the factor 2 uncertainty in the absolute space den- 
sity normalization of our distributions. The dashe d curve shows the b ivariate 
distribution function of the form presented in de Jons & Lacev 1 200C ) fitted 
to the data. The dashed histogram is the space density distribution derived 
bv lBlanton et al] 12003) from SDSS data for galaxies with n > 3 and scaled 
as described in the text. 



Figure fT2l also shows the fitted bivariate distribution function 
of the form described by Ide Jo ng & L acevl l2000l) . This function 
takes the shape of a lSchechted \ 19761) luminosity function in the 
luminosity dimension (with the usual parameters (/>,, M*, and a), 
and of a log-normal distribution in the scale size dimension, with 
a width a\, a median scale size for galaxy with M = M* of R et , 
and a shift in median scale size with luminosity parametrized by j3: 



<f>(M, log(Re))dMdlagRe = 0.41n(10) 



ln(10) 



27T<7 A 



(6) 



;^ 1q -o.4(«+i)cm-m.) ex p(-10-°- 4(M - M * ) )dM 



x exp 



■*( 



log(fl c /fl e »)-0.4j3(M-A/,) 
<r A /ln(10) 



d log R c 



The second line in this equation represents the Schechter luminos- 
ity function in magnitudes, while the third line represents the log- 
normal scale size distribution. 

The function was fitted using a maximum likelihood tech- 
nique and the resulting parameters are listed in TableQ with the 
errors indicating the 95% confidence limits. These confidence lim- 
its were determined by means of Monte Carlo bootstrap resampling 



of the galaxies IPress et aT]|l993h . where we also varied the size 
of the bins in both directions. The fitted function is a good de- 
scriptio n of the o bserved distribution: the goodness-of-fit param- 
eter Q (Press et al. 1993) was larger than 0.1 in more than 28% of 
the bootstrap resampled realizations and Q > 0.001 in more than 
87% of the realizatio ns. A function similar in shape was fitted by 
ICholoniewskil ^)85) to a sample of 233 elliptical galaxies and our 
parameters do agree to within the quoted uncertainties. 

In Fig. | 12 | we compare our bivariate distribution with the find- 
ings of Blanton et al. (2003), based on the SDSS data (see their 
Fig. 13). They fitted Sersic profiles to the total galaxy luminosity 
profiles, without an added exponential disc component and hence 
without doing a full bulge/disc decomposition. We show their 01 i 
distribution for all galaxies with an Sersic n > 3. We scaled their 
distribution to our Hu bble constant and us ing a 0,1 i to R conver- 
sion of 0.2 mag (from lBlanton et alj|2003l 01 (r - i) 



2.51og(l.l), and R 



0.3 from lBell & Vandenberd 



-0.4, 



1198 

Surprisingly, the distributions agree well at the bright end, 
even though our absolute space density zero-point calibration is 
somewhat uncertain (see errorbar in top panel) due to the clus- 
ter over density in the EFAR region and our limited known- 
ledge of cluster inco mpleteness. The cutoff at ii o >20kpc in the 
iBlanton et alj 120031) distribution is not real, but due to the cut- 
off in their figure. At fainter luminosities our distribution is lower 
than that of IBlanton et alJ 1200 3h . One explanation for this might 
be that we were stricter in excluding spiral galaxies from our dis- 
tribution, while the Blanton et al. distribution includes all galaxies 
with n > 3, which may include spiral galaxies whose combined 
bulge and disc luminosity profile mimics that of an elliptical galaxy. 
On the other hand, incompleteness may affect our distribution at 
the faint end. Our sa mple was drawn from rich clusters, while the 
IBlanton et al]l2003t) SDSS distributions were calculated for all sky. 
This will have little effect on the bright end of the distribution, be- 
cause the brightest elliptical galaxies are only found in cluster, but 
will affect the faint end. Notwithstanding, the size distributions at 
a given luminosity are very similar. 

In recent yea rs it has become clear there are two types o f 
elliptical galaxies iKormendv & B ender ll996tlFaberetalJll997t) . 
On the one hand we have the bright elliptical galaxies with boxy 
isophotes, with a clear central core, dynamically mainly pressure 
supported, while on the other hand we have the lower luminosity el- 
liptical galaxies, often with more discy isophotes, cuspy cores and 
for a larger fraction supported and flattened by rotation. The boxy 
galaxies tend to be central dominant galaxies of clusters (cD galax- 
ies) or at the center of one of the subclumps of a cluster (resulting 
from cluster merging?). The transition from one type into the other 
occurs at about -21.2 l/-mag iFabe r et alj[l997h . i.e. about -20.5 
i?-mag, with an overlap of both types occurring for about 1.5 mag. 

The bright, boxy type elliptical galaxies are probably the re- 
sult of the merging of many subclumb units in the hierarchical 
galaxy formation scenario. The origin of the lower luminosity el- 
liptical galaxies is less clear. They could be t he result of th e merg- 
ing of two similar sized gaseous galaxies (Barnes 1992) or the 
result of "failed" spiral galaxies due to the lack of angular mo- 
mentum in the proto-galaxy ( e.g., IBalcanton. Spergel. & Summers 
Il997t lde Jong & Lacevl 2000). These two options might actually be 
somewhat related, as merging of two sub units is of course the result 
of lack of orbital angular momentum. For SO galaxies even more 
formati on scenarios have been proposed: ram-pressure str ipping 
of ga s <Gunn & Got3ll972t lKentlll98lt lOuilis. Moore. & B ower 
2000), galaxy harassment via high speed impulsive encounters 
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Table 1. Bivariate distribution function parameters. Errors indicate 95% confidence limits. 



Fit 


</»* 


a 


M* 




0"A 


P 




(xfCT 4 Mpc- 3 ) 




(i?-mag) 


(kpc) 






Total Galaxy 


0.99 ±0.19 


1.06 ±0.26 


-21.81 ±0.18 


3.9 ±0.5 


0.34 ± 0.04 


-0.62 ± 0.08 


Bulge 


0.83 ±0.31 


1.30 ±0.50 


-21.15 ±0.30 


2.3 ±0.4 


0.39 ± 0.04 


-0.55 ±0.10 


Disc 


0.91 ±0.26 


-0.17 ± 0.24 


-21.98 ± 0.35 


12.6 ±2.2 


0.42 ± 0.04 


-0.40 ± 0.07 
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Figure 13. The space density of EFAR galaxies as function of magnitude 
and effective radius of the total galaxy (solid histogram). Errorbars and up- 
per limits the same as in Fig. 1 121 The dashed histogram is the space density 
distribution derived by Blanton et al. 1 2003) from SDSS data for galaxies 
with n < 1 . 5 and scaled as described in the text. 



iMoore et al]ll99fJ) . cluster tidal forces jBvrd & Valtonenlll99fi) 
which distort galaxies as the y come close to the c entre, in- 
teraction/merging of galaxies llckdll985l iBekkil Il998h . and re- 
moval & co nsumptio n of the gas due to the cluster environ- 
ment lLarson. Tinslev. & Caldwel3ll980l : lBekld. Couch. & Shioval 
120021) . In all these scenarios the origin of non-boxy elliptical galax- 
ies and SO galaxies lies in spiral galaxies or proto-spiral clumps. 
Therefore, it seems natural to compare the size-luminosity distri- 
bution of our ellipical sample to that of spiral galaxies. 

In Fig.[21 we show again the space density distribution of 
EFAR galaxies as function of size and luminosity, bu t this t ime we 
compare the distribution with that of Blan ton et alJ {2003) SDSS 
galaxies with Sersic profile fits with n < 1.5, i.e. galaxies with 
exponential-like profiles, predominanlty spiral galaxies. The Blan- 
ton et al. distributions were scaled for Hubble constant and pass- 
band in the same way as in the previous diagram for n > 3 galaxies. 
The disc effective radii were multiplied by 0.75 to take into account 
that for an R * -law galaxy the projected surface brightness effective 
radius is about 0.75 the spherical effective radius. In this way, we 



make a rough comparison between the stellar mass distributions, 
ignoring the effect of other baryonic components (most notably the 
HI in the disc galaxies, which could be converted to stars during 
merging) and the systematic change in M/L with luminosity, es- 
pecially for disc galaxies. When scaled for Hubble const ant and 
passba nd, the Sc-Sd spiral galaxy d istribution of | d e long & Lacevl 
(2000) is very similar to the one of iBlanton et all jlfJOsT Tparticu- 
larly at the bright end. 

Figure lT3l immediatelv shows that to form the brightest ellipti- 
cal galaxies one will need to merge several current day bright spiral 
galaxies, because there is a substantial space density of bright el- 
liptical galaxies that are at least 0.75 magnitude brighter than the 
brightest spiral galaxies. In the process the effective radius of the 
resulting galaxy is increased by a factor of few compared to the 
originating galaxy, presumably because the kinematic energy of 
the galaxies orbiting each other is converted to h eat the r andom 
kinematic energy of the stars in the merger product, barnesl il992h 
finds for dissipationless equal mass mergers a typical increase of 
effective radius of about 1.4, similar to what on e expects assum- 
ing conservation of energy and the virial fheorum IColeet all 2000: 
IShen et aljEooi) . We have stepped our luminosity bins in Fig. l 131 
by 0.75 mag, i.e. about a factor 2 in luminosity (and roughly stel- 
lar mass), so each equal mass dissipationless merger of two disc 
galaxies would result in a merger remnant one bin up with about 
1.4 larger size. Merging of at least a few and probably many more 
of present day spiral galaxies would be needed to form the brightest 
elliptical galaxies. Because merging of spiral galaxies will involve 
dissipation and therefore the scale size may be less, it may be more 
reasonable to think that the brightest elliptical galaxies formed from 
the merging of several smaller E/S0 galaxies. 

A similar conclusion was reached bv lShen et alJ l2003). They 
showed that by having a series of mergers of small elliptical galax- 
ies they could create the observed size-luminosity relation and its 
dispersion, when making reasonable assumptions about the trans- 
fer of orbital angular energy. However, they did not compare space 
densities, and furthermore, this process leaves the question of the 
origin of small elliptical galaxies unanswered. 

The distribution of the lower luminosity elliptical galaxies 
(A/ to t>-22.5 i?-mag) in Fig.^|is more surprising. Even though in 
number density not completely unlike the spiral galaxies (remem- 
ber the absolute zero-point for the elliptical distribution is some- 
what uncertain), the effective radius distribution of the elliptical 
galaxies peaks at smaller radii than the distribution of disc galax- 
ies. As argued before, dissipationless merging results in larger scale 
sizes, typically by a factor of about 1.4, and therefore the dissipa- 
tionless merging of typical current day spiral galaxies cannot result 
in these current day lower luminosity elliptical galaxies. 

What other creation scenarios can we invoke to explain the 
lower luminosity elliptical galaxies? First of all, the merging of two 
spiral galaxies is not going to be dissipationless. Simulations that 
include gas show that during the merger process gas will quickly 
stream to the center of the merging galaxies, presumably creating a 
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starburst and dragging along some of the existing stellar population, 
creating a more concentrated stellar remnant. Whether this process 
is enough to offset the effect of the added orbital kinematic energy 
is questionable, given that current bright spiral galaxies have gas 
mass fractions of 5-20%. At higher redshifts the gas fractions will 
have been higher and therefore this process may have been more 
efficient then. In the hierarchical galaxy formation picture, galaxies 
at higher redshifts have smaller scale sizes, so the small scale sizes 
of low luminosity elliptical galaxies in clusters may be the result of 
high redshift mergers which have subsequetially not been able to 
grow any more in the cluster environment. 

Another option to explain the small scale sizes of lower 
luminosity elliptical galaxies might come from the failed disc 
scenario propos ed by several authors fo r different reasons 
(e.g. iDalcanton. Spergel. & Summers! Il997t iMo. Mao. & White! 



19981 Ide Jong & Lacevl l200Ct iMcGaugh & de Blokl Il998l) 
de Jong & Lacev 1 2000) found that the scale size distribution of 
spiral galaxies used here was narrower than what one would expect 
in hierarchical galaxy formation scenario where disc galaxy sizes 
are determined by their spin acquire d from tidal torques with 
neighbors. They proposed (similar to IMo. Mao. & White! f 1998) 
that the lowest angular momentum proto galaxies never were able 
to form a proper disc and immediately collapsed to spheroids. 
The problem here is that we cannot compare the elliptical galaxy 
distribution drawn from a cluster sample to the distribution of 
spiral galaxies drawn mainly from the field. 

Can spiral galaxies falling into a cluster result in a elliptical 
scale size distribution as observed here? Simple stripping of gas by 
ram-pressure or tidal interactions and subsequent fading of the stel- 
lar population is clearly not enough, elliptical galaxies are just too 
small compared to spiral galaxies. Subsequent stripping of the stars 
in the o uter parts of the infalling disc galaxy by galaxy harassment 
(Moore etallll996l) could result in smaller effective radii of the 
remnants. Still a lot of stripping would have to occur, as stripping 
half the stellar mass exclusively from the outer part of the galaxy 
would shift the disc distribution one bin down while dividing the 
size by two. A considerable fraction of this process would have to 
occur at quite high redshifts (z > 1.5) as only then bright disc 
galaxies were substantially smaller th an current day disc ga laxies 
Ide Jo ng & Lacevlll99 9: Ravindranath et al. 2004; Fergus on et alj 

As so often, several of the options listed above combined may 
be the final answer. However, these relative elliptical and spiral 
galaxy scale size distributions suggest that a substantial fraction of 
the lower luminosity elliptical galaxies were created at quite high 
redshift, not only in stars but also in structural parameters, when 
disc galaxies were still smaller and more gas rich. 

Until so far, we have only looked at the luminosities and sizes 
for the total systems, but in Figs ll4l and ll5l we show the luminosity- 
size distributions for the bulges and the discs, respectively. The lu- 
minosity and size parameters are much better determined for the 
total system than for its components, but nontheless it will be in- 
structive to compare the distributions of the components to that of 
spiral galaxies. 

Comparison of the elliptical bulge parameters with the SDSS 
disc parameters in Fig. 1141 shows that the distributions are some- 
what similar, but that the bulges are about a factor ~1.5 to 2.5 
smaller than the spiral discs from the brightest to the faintest lu- 
minosities. The scaling of the SDSS disc parameters is here the 
same as in Fig. 1131 including the factor 0.75 decrease in disc scale 
sizes to account for the difference of a spherical bulge projected to 
a 2D effective radius versus a 2D disc effective radius. 
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Figure 14. The space density of EFAR elliptical galaxies as function of 
bulge magnitude and bulge scale length (solid histogram). Errorbars and 
upper limits are the same as in Fig. 1121 The dashed histogram is the space 
density distribution derived bv lBlanton et al] 1200 3l) from SDSS data for 
galaxies with n < 1.5, scaled as described in the text. 



On the other hand, when we compare in Fig.HU the SDSS 
n < 1.5 disc parameter distribution with the distribution for the 
component identified as discs in our EFAR sample, we see again 
somewhat similar, parallel distributions, but this time the EFAR 
disc-like component is about a factor two larger than the SDSS 
discs. We have to be careful, because as argued before in section l33l 
when discussing Fig.0 discs with low surface brightness and/or 
small scale length might be hard to detect when superposed on a 
high surface brightness bulge and therefore this distribution may 
be incomplete. Still, it is clear that there are larger disc components 
in the EFAR non-spiral galaxies than seen in the SDDS spiral discs. 
Upper limits are somewhat hard to draw on this diagram, as the se- 
lection was dominated by the bulge light distribution, not by the 
disc light. 

Without kinematic information proving that the B/D decom- 
positions we made of our sample has an underlying physical mean- 
ing it is somewhat premature to interpret these distributions in de- 
tail, but we can make some speculations. Could the spheroidal com- 
ponents be the result of early, disipationless mergers of roughly 
equally sized galaxies from the time when galaxies were still much 
smaller (z>l-2)? The disc-like component can clearly not be 
formed from normal disc galaxies, there scale sizes are just too 
large for their luminosity. The disc-like component tends to dom- 
inate a bit in the most luminous galaxies (see Fig.[5). Are these 
the stars from smaller, infalling galaxies which maintain a larger 
component from their initial orbital energy? Do they form disc-like 
components because they predominantly fall into the cluster on the 
central galaxies along the filamentary structures of the cosmic web? 

While attractive from the galaxy structure point of 
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Figure 15. The space density of EFAR galaxies as function of disc mag- 
nitude and disc scale length. Errorbars the same as in Fig. 1121 The dashed 
histogram is the space density distribution derived by Blant on et al]l2003t) 
from SDSS data for galaxies with n < 1.5, scaled as described in the text. 
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Figure 16. The space density as function of bulge magnitude and bulge 
effective radius for the EFAR galaxies (solid histogram) and for the semi- 
analytic models of Cole et al. 12000) (dotted histogram). Errorbars and up- 
per limits the same as in Fig. 1121 



view, this picture leaves some other properties of ellipti- 
cal galaxies unexplained. Elliptical galaxies in clusters obey 
a colour-magnitude relation up to r edshifts of ~0.9 (e.g., 
Stanford. Eisenhardt. & Dickinson! 1 1 9981) . which means that this 
relation is driven by a mass-metallicity relation, bigger galaxies 
being more metal rich. This therefore means that the proto-galaxies 
that formed the big spheroids in the above picture would have to 
have higher metallicity before they merged than the proto-galaxies 
of smaller galaxies. Big elliptical galaxies also have higher a- 
element over Fe abundance ratios than smaller galaxies, meaning 
that th eir star forma tion time-scales must have been shorter (e.g., 
iMehlert et ai]l2003l) . In the hierarchical galaxy formation picture 
this could be accomplished by having the proto-galaxies that will 
merge into big galaxies form early on top of an over-density of 
what will later become a cluster of galaxies at redshifts much larger 
than 3. In a ACDM Universe there would be a lot of early merging 
causing vigorous and short time-scale star formation, leading to 
non-Solar rates of a-element over abundances. Because these 
proto-galaxies would already live in the deeper potential well of 
what will later become a large galaxy in a cluster of galaxies, they 
would be better in holding on to their enriched gas and therefore be 
able to create the mass-metallicity relation. However, in the above 
picture galaxies falling in along the filaments would have lower 
(over-)abundances and would potentially be somewhat younger, 
and hence t he disc components should have lower metallicity and 
age. While Ide Jong & Daviesl J 19971) show that discy elliptical 
galaxies indeed have line indices indicating younger ages and/or 
lower metallicities, spatial resolved line indices mapping as per- 



formed by for instance the SAURON team iFalcon-Barroso et alj 
2004) will be needed to fully access this model. 

Figure fTol shows the comparison of the (bulge magnitude, 
bulge effective radius)-bivariate spac e den s ity of the EFAR ellip- 
tical galaxies with the models o f lCole et all (2000). The models are 
over-predicting the effective radii of early-type bulges by about a 
factor of two, even more at fainter magnitudes. This is the same 
problem as identified in Fig. ll 31 the effective radii of low luminos- 
ity early-type galaxies are smaller than those of spiral galaxies of 
similar luminosity. Cole et al. assumed that all early-type galax- 
ies were created by merging in their models, and as their models 
do a reasonable job in modelling the spiral galaxy scale size dis- 
tribution Ide Jong & Lacevl2 000). the models will over-predict the 
early-type galaxy scale size in the same way as was discussed for 
HgJSI 



4 SUMMARY AND CONCLUSIONS 

We have performed 2D bulge/disc decompositions on a sample of 
558 early-type galaxies from the EFAR sample using the GIM2D 
package. In contrast with most earlier work, we have used a Sersic 
luminosity profile for the spheroidal component, while using the 
common exponential light profiles for the disc-like component. 
We showed with extensive testing on model galaxies that the to- 
tal galaxy parameters and disc parameters can be recoverd at the 
10-15% level, but that the bulge parameters are only recovered at 
the 30% level, mainly due to degeneracies resulting in correlated 
errors between ub and R c . Our galaxy selection criteria are well 
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known and we have calculated bivariate space densities of elliptical 
galaxies as function of galaxy parameters. 

Without kinematic or population information it is impossible 
to prove that our bulge and disc parameter are physically meaning- 
ful, but under the assumption that the spheroid and disc parameters 
represent real galaxy components we derive the following conclu- 
sions: 

(i) The scale sizes of the bulge and disc components are corre- 
lated. This may be partly due to the exclusion of spiral galaxies 
from the sample, excluding objects with small bulges and large 
discs. Howev er, it should be noted that sp iral gala xies show a 
similar trend Ide Jondfl 996b: [Cou rteau. de Jong. & Broeilslll996t 
Mac Arthur. Cour teau. & Holtz marl2003l) . 

(ii) The bulge-to-disc B /T ratios, the Sersic ub bulge shape 
parameters, and the bulge effective radii show positive correlations 
with each other. We have no model that tells us which is cause and 
what is effect. 

(iii) The median ub value is 3.24 for all galaxies in our sample, 
is 3.66 for all non-spiral, bulge dominated galaxies in our sample, 
and is 4 for non-spiral galaxies with B/T>0.1. Given that most 
spiral galaxies have bulge ns values much lower than these values, 
this means that the standard Ba de Vaucouleurs (i.e. n_e=4) fitted 
to spheroids may not be the most appropriate choice when doing 
bulge/disc decompositions on a random set of galaxies. In cases 
where one does not have the signal-to-noise or the spatial resolu- 
tion to perform full Sersic bulge plus exponential disc fit, it will be 
more correct to use fixed value of ns~3.5 for the brightest galax- 
ies. Lower values should be used for less luminous galaxies. 

(iv) The (luminosity, scale size) bivariate space density distribu- 
tion of bright elliptical gala xies is well described by the analytic 
parametrization presented bv lde Jong & Lacevl 42000). 

(v) Comparing the total galaxy (luminosity, scale size) bivari- 
ate distributions of EFAR non-spiral galaxies with SDSS disc-like 
galaxies shows that the brightest early-type galaxies could in prin- 
ciple be formed by merging a few large, current day spiral galax- 
ies. Low luminosity early-type galaxies have much smaller scale 
sizes than spiral galaxies of the same luminosity, and hence they 
can surely not be created from simple merging of two current day 
spiral galaxies. 

(vi) When comparing the EFAR bulge and disc (luminosity, 
scale size) bivariate distributions to those of the SDSS disc-like 
galaxies, we showed that the bulges of early-type galaxies are typ- 
ically a factor 1.5 to 2.5 smaller and the discs a factor 2 larger than 
current day disc galaxies at a given luminosity. 

We speculate that the spheroidal components are the result of 
merging of similar sized small proto-galaxies galaxies at high red- 
shifts, while the disc-like component may be the result of smaller 
galaxies falling in later along the filaments of the cosmic web. 
While attractive from the galaxy structure point of view, this model 
requires careful tuning of the formation process of the proto- 
galaxies at high redshift in order to reproduce the mass-metallicity 
relation and the a-element over-abundace observed at lower red- 
shifts. 
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